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Abstract 

In this paper, a recently developed numerical method to analyze dielectric-spectroscopy 
data is applied to a-phase polyvinylidene fluoride (PVDF). The numerical proce- 
dure is non-parametric and does not contain any of the extensively used empirical 
formulas mentioned in the literature. The method basically recovers the unknown 
distribution of relaxation times of the generalized dielectric function representa- 
tion by simultaneous application of the Monte Carlo integration method and of the 
constrained least-squares optimization. The relaxation map constructed after the 
numerical analysis is compared to a-phase PVDF data presented in the literature 
and results of the parametric analysis with a well-known empirical formula. 
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integration, constrained least-squares. 



Dielectric spectroscopy is a well established tool for materials characteriza- 
tion With this technique, the complex dielectric permittivity ^(a;) is ob- 
tained as a function of angular frequency u {u = 2txv] v is the frequency of 
the applied voltage). Analysis of ^(ci;) is usually performed with empirical iox- 
mulas 0, 0, 0, 1^ . A neglected approach is the distribution of relaxation times 
(DRT), which is based on functional analysis (it requires the solution of an 
integral equation) i. The DRT expression is a Volterra equation [2|, which 
is a special form of the Fredholm integral equations [1]. Such equations are 
usually considered to be ill- conditioned because of their non-unique solutions. 
However, the approach used here and recently presented elsewhere jl, 0,113, U 



leads to unique solutions by means of a constrained least-squares fit and the 
Monte Carlo integration method. Macdonald [121 stated that approaches with 
summations of delta functions would not suffer from the limitations of ill-posed 
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Fig. 1. (a) Real and (b) imaginary parts 
of a-phase PVDF dielectric function of 
temperature at various frequencies v: (o) 
0.1 Hz, (□) 1 Hz, (>) 10 Hz, (^) 100 Hz, 
(<) 1 kHz, (*) 10 kHz, (v) 100 kHz and 
(A) 1 MHz. The lines are drawn to guide 
the eyes. 
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inversions if complex nonlinear least-squares are employed, another approach 
to solve inverse problems. In this paper, we apply this numerical procedure 
that was developed by one of us (ET) to a semi -cry stalhne polymer, a-phase 



polyvinylidene fluoride [a- phase PVDF) [liQIilES- 

An empirical formula 

and results from the literature are also included in order to evaluate the new 
DRT approach and to compare several analyzing techniques. Impedance mea- 
surements on a-phase PVDF are performed at various tempatures T between 
173 K and 450 K within frequency windows of 0.1 Hz < z/ < 1 MHz. 

In Fig. 1, the relaxation spectra of a-phase PVDF are shown at several fre- 
quencies as a function of temperature. Local movements of the structural 
units in a-phase PVDF are observed at low temperatures, T < 250 K. This 
relaxation is usually called /3 relaxation. In a-phase PVDF, the magnitude 
of the /3-relaxation is much smaller than those of two a-relaxations. The a- 
relaxations are related to cooperative movements of segments of the macro- 
molecular chains. The subscripts 'a' and 'c' denote the fast asymmetric (amor- 
phous) and the slow symmetric (crystalline) relaxation modes in the mate- 
rial |l3|, |l^ , respectively. 



The dielectric function of a material can be expressed in a general forma as 
follows, 



£(^) = e + x(^) + c^/(^£o^), (1) 

where e and a are the frequency-independent relative dielectric permittivity 
and the ohmic conductivity of the material, i = and eo = 8.854 pFm~^. 
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x(i^) is the susceptibility of the material which is usually expressed through 
the Havriliak-Negami empirical formula (HNEF) as 



xM = EA^.[i + (^^^.rr"^- (2) 

Here, we consider that there are several relaxation processes in the mate- 
rial (summation over j). The parameters Tj and Ae^ are the most probable 
relaxation time and the magnitude of the related polarization relaxation, re- 
spectively. The parameters aj and /5j determine the shape of the relaxation, 
aj < 1 and ajPj < 1 [l^. A curve- fitting procedure based on the complex- 
nonlinear-least-squares (cnlsq) method would yield the four unknowns for 
each relaxation {r, Ae, a, /?}, as well as the high-frequency dielectric permit- 
tivity e and the ohmic losses proportional to a 

In the DRT representation, the susceptibility x{^) of Eq. (1) is expressed as 



x{to) = Ae g(r)(l+z^r)-Mr, (3) 

JO 

where g(r) is the distribution of relaxation times. Eq. (3) is the Volterra 
equation |7i and (1 + iujt)~^ its kernel )C{uj, r). For the numerical procedure, 
several relaxation times are selected from a log-linear scale on the frequency 
axis 1^ . The number k of the relaxations is smaller than the total number 
n of experimental data points. In that case, the integral equation in Eq. (3) 
becomes a sum that is linear in the gfc. 



k 



(4) 



Since the are expected to be positive, gfc > 0, a constrained least-squares 
algorithm is employed with newly selected r sets for N times, a-phase 

PVDF had previously been analyzed by Grimau et al. 14| and by Bella et al. 
[lil . They used a Metropolis algorithm to obtain the relaxation characteristics 
of a-phase PVDF. 



Below, we first present and compare the results obtained from HNEF and DRT 
analyses at 173 K. In Fig. 2, HNEF analysis is plotted. The symbols in Fig. 2a 
are the experimental values, and the line is the curve fitted with Eq. (2) 
for j = 2. The relative error 6 of the real and imaginary parts of e{u!) are 
presented in Fig. 2b. We illustrate the relative error in order to compare the 
DRT analysis with other methods. The average magnitude of the error in the 
imaginary part 6e" is approximately ten times higher than that of the real part 
6e', 6e" ~ 10 6e'. In addition, a detailed analysis of the data shows that 90% 
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Fig. 2. (a) Cole-Cole parametric plot of the dielectric function at 173 K. The symbols 
and the line represent the experimental data and the curve obtained from cnlsq 
with the Havriliak-Negami empirical formula, respectively, (b) Relative error 5 of 
the CNLSQ fitting. 

and 80% of the data fall into the region Se" < 58 x 10"^ and Se" < 48 x 10"^, 
respectively. 

In Fig. 3, the DRT analysis based on Eq. (3) is presented. Here, the data are 
divided into three regions (the results are presented in Fig. 3a and 3b with open 
symbols). The number of relaxation times in each minimization step is chosen 
to be 2/3 of the number of the experimental data points. The Monte Carlo 
loop A^, which is the number of minimizations performed, is 2^^. The average 
magnitude of the error in e" is again approximately ten times that of the 
real part, Se" ^ lOSe', as shown in Fig. 3b. Compared to the HNEF, the DRT 
analysis yields lower relative errors. 90% and 80% of all the data points fall into 
the region 6e" < 40 x 10~^ and Se" < 26 x 10~^, respectively. These values 
are nearly 50% lower than the corresponding values of the HNEF analysis. 
In addition, the distribution of the errors is more symmetric than with the 
HNEF analysis. The extracted DRT spectrum is presented in Fig. 3c, which 
is plotted with error bars after presenting the data in bins. The numerical 
procedure predicts smooth continuations of the relaxation times r outside the 
experimental window. The experimental frequency window is indicated with 
dashed vertical lines (- - -) in Fig. 3c. 



Application of the DRT analysis to data at different temperatures yields dis- 
tributions similar to that presented in Fig. 3c. Levy statistics j2l[ 2^ | will 



be later used to assess these distributions (see Refs. [10, lUl). The mean re- 
laxation times from the Levy statistics and from the HNEF analysis are em- 
ployed to calculate the temperature dependence of the relaxations with the 
Vogel-Fulcher-Tammann (VFT) expression 



r = Toexp{W[h{T + To)]-'} (5) 
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Fig. 3. (a) Cole-Cole parametric plot of the dielectric function at 173 K, the symbols 
and the line represent the experimental data and the curve obtained from the DRT 
analysis, respectively. The frequency window is divided into three regions on the 
frequency- axis, which are indicated with different symbols ('o';high, '□';intermediate 
and '[>';low). (b) Relative error 5 of the DRT analysis, (c) The DRT of a-phase 
PVDF for the experimental window 0.1 Hz < < 1 MHz; the dashed vertical lines 
(- - -) indicate the lower and upper limits of the relaxation times r probed in the 
experiments. 

where tq is a constant, W is a constant related to activation energy of the 
relaxation process, and Tq is the Vogel temperature related to the glass tran- 
sition Tg of the material. For Tq = 0, Eq. (5) reproduces Arrhenius' expression, 
r = To exp[W {kbT)~^] with kb = 86.1321 |j,eVK~^ In Fig. 4, the relaxation map 
extracted from the DRT analyses at different temperatures is shown with gray 
scale. In order to compare both analyses, we plot the DRT and HNEF results 
with solid ( — —) and dashed (- - -) lines, respectively. Moreover, results from 
the literature 3, lil are also given [chained line (- ■ -)]. The (3 relaxation is 
omitted, since it has a very broad distribution and a low magnitude compared 
to the a-relaxations as presented in Fig. 3c. The fit parameters for the VFT 
and Arrhenius expressions and the values from the literature are presented in 
Table 1. The relaxation is known to be asymmetric, and our investigations 
by means of Levy statistics confirm this observation — we are able to fit two 
symmetric relaxations that result in two different VFT expressions, see the 
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Table 1 

Fitting parameters for the a relaxations. Tg is determined by the temperature at 
which the VFT fits yield logr = 2. 





log (to [s]) 


W [eV] 


To [K] 


Tg [K] 


Note 


tta relaxation 


This work 


-13.6 


0.079 


191.0 


223.4 


DRT 


This work 


-12.8 


0.085 


191.3 


224.0 


DRT 


This work 


-13.1 


0.076 


195.6 


225.4 


HNEF 


Bella et al. [1^ 


-13.1 


0.088 


187.5 


222.1 


HNEF 


Oc relaxation 


This work 


-16.6 


0.40 






DRT 


This work 


-18.4 


0.45 






HNEF 


Bella et al. \13] 


-19.3 


0.49 






HNEF 



solid lines ( ) in Fig. 4. The distribution is symmetric at low temper- 
atures, but becomes asymmetric around 250 K. When the fitting results are 
taken into consideration, there are slight differences between our studies and 
those of Bella et al. [l^ and Grimau et al. [3], which are probably due to 
differences in the preparation and the thermal history of the materials. 



Here, we presented the results of dielectric-data analyses for a-phase PVDF 
with two different techniques. For our particular material, both techniques 
result in equivalent temperature dependencies, which were also compared to 
data in the literature. Our numerical method for finding the DRT yields unique 
distributions for given dielectric data without any a-priori assumption. We 
have thus demonstrated that the dielectric data of materials can be analyzed 
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Fig. 4. Relaxation map for a-phase 
PVDF, which clearly shows the aa and 

Oc modes. Solid ( ) and dashed 

(- - -) lines represent the fits to VFT or 
Arrhenius expressions for the DRT and 
HNEF analyses, respectively. The results 
of Grimau et al. ^IJ,] and Bella et al. 
[l^ are plotted as chained lines (- • -). 
The gray scale illustrates the probabil- 
ity density of the relaxations g(T), dark 
areas are the most probable regions for 
the particular relaxations. The omitted 
(white) region, labeled R, would contain 
molecular information related to melt- 
ing of the sample. 
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with a more general representation than the empirical formulas. The DRT 
method might be useful to investigate the origin of dielectric relaxations and 
the electrical properties of materials. 
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